Many-body exchange-correlation effects in graphene 
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We calculate, within the leading-order dynamical-screening approximation, the electron self- 
energy and spectral function at zero temperature for extrinsic (or gated/doped) graphene. We 
also calculate hot carrier inelastic scattering due to electron-electron interactions in graphene. We 
obtain the inelastic quasiparticle lifetimes and associated mean free paths from the calculated self- 
energy. The linear dispersion and chiral property of graphene gives energy dependent lifetimes that 
are qualitatively different from those of parabolic-band semiconductors. 
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INTRODUCTION 

Recent developments in techniques for fabricating con- 
ducting graphene layers [l[ have provided the physics 
community with a unique opportunity to study an in- 
teracting two-dimensional (2D) massless Dirac fermion 
system. This has led to considerable experimental and 
theoretical activity in this field The band structure 

of graphene, by dint of its honeycomb lattice, has linear 
dispersions near the K and K' points of the Brillouin zone 
Q . The corresponding kinetic energy of graphene for 2D 
wave vector k is given by e^ s = S7|k|, where s = ±1 indi- 
cate the conduction (+1) and valence (—1) bands, respec- 
tively, and 7 is a band parameter (and Fermi velocity of 
graphene is given by vp = "f/h). The corresponding den- 
sity of states (DOS) is given by D{e) = g s g v \e\/ {2ir"f 2 ) 1 
where g s = 2, g v = 2 are the spin and valley degen- 
eracies, respectively. The Fermi momentum (kp) and 
the Fermi energy (Ep) of 2D graphene arc given by 
kp = (Ann) g s g v ) x l 2 and Ep = jkp where n is the 2D 
carrier (electron or hole) density. 

Many electronic properties of a system are strongly in- 
fluenced by the presence of electron-electron interaction. 
In this paper, we investigate theoretically the electron- 
electron interaction induced exchange-correlation effects 
in two dimensional graphene layer. We calculate, within 
the leading-order dynamical-screening approximation, 
the electron self-energy for extrinsic graphene at zero 
temperature. The self energy is the central quantity that 
determines the other Fermi liquid parameters. We ob- 
tain the single particle spectral function and the inelas- 
tic quasiparticle lifetimes and associated mean free paths 
from the calculated electron self-energy. 

The self-energy is given by in the screened interaction 
(GqW approximation) [5| 
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X/ksT, s,s' = ±1 denote the band indices, 
) = l/(iu/ n — £ks) is the unperturbed Green's 



function (£k s = e^ s — fi where /i is the chemical poten- 
tial), V c (q) — 2ite 2 / nq is the bare Coulomb potential 
(with background dielectric constant k), e(q,u>) is the 
dynamical screening function (dielectric function) given 
by e(q, w) = 1 — v c (q)H(q,u), where U(q,uj), the 2D po- 
larizability, is given by the bare bubble diagram [fj. In 
Eq. ([T]) F ss / (k, k + q) is the overlap of states, given by 
F ss r(\t, k') = (1 + ss' cos #kk')/2, where 9 is the angle 
between k and k'. 

After the standard procedure of analytical continua- 
tion, the self-energy can be separated into the exchange 
and correlation parts £ a (k,£j) = Sf(k) + S™ r (k, w). The 
exchange part is given by 
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where n_F(£ks) = 0(£ks) is the Fermi function at T = 0. 
The correlation part, X^ or (k, oj), is defined to be the part 
of E s (k, ui) not included in £° x . Since the exchange part 
of self-energy E° x (k, w) is studied in ref. Q we provide 
in this paper the results of correlation part. In the GW 
approximation, the £™ r (k, to) can be written in the line 
and pole decomposition S™ 1 ' = E|, inc + £P° lc , where 
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The Sii ne is completely real because e(q, iu) is real. Thus, 
ImfEpoie] gives the total contribution to the imaginary 
part of the self-energy. In this calculation we use the 
cut-off at the wave vector k c with respect to the Dirac 
point; i.e., the integral is cut off at k — k c . From the 
calculated self energy H(k, u>), we can obtain the single- 
particle spectral function given by 
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FIG. 1: (a) Real (top panel) and imaginary (bottom panel) 
parts of the self-energy u) and (b) spectral function 
A(k,u>) as a function of the energy uj for k = (the Dirac 
point). 
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FIG. 2: (a) Real (top panel) and imaginary (bottom panel) 
parts of the self-energy E(fc, uf) and (b) spectral function 
A(k,ui) as a function of the energy ui for k — kF (at the 
Fermi surface). The vertical arrow in (b) denote 8- function 
in the spectral function of weight (2tt x 0.9) at uj = 0. 
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The spectral function A(k,u) satisfies the sum rule 

/OO J 
_A(fc,o;) = l J (6) 

which is generally satisfied to within less than a percent 
in our numerical calculations. 

RESULTS 

In Figs. 1 and 2 we show the real and the imaginary 
parts of the self-energy and spectral function for k = 
(the Dirac point) and k = Uf (Fermi surface). Through- 
out this paper we use the following parameters: dielectric 
constant K — 2.5, Fermi velocity vp = j/H = 10 8 cm/s, 
and cut-off momentum k c = l/a, where a = 2.46 A is 
the lattice constant of graphene. 

For k = (Fig. 1) we find a strong peak in Im[S(fc, w)J 
(associated with a finite step in Re[£(/c, lj)]) due to plas- 



mon emission. Since Im[S] is finite for all energies except 
lu = we expect a finite damping in the spectral func- 
tion. In spectral function we find two peaks. The first 
peak (from lu = 0) in A(k,u>) corresponds to the usual 
quasiparticle (i.e. a bare particle surrounded by a cloud 
of virtual plasmons and particle-hole excitations), and 
the second peak corresponds to a plasmaron, which is in- 
terpreted as a hole coupled to a cloud of real plasmons 

1- 

In Fig. 2 we show the self-energy and spectral function 
for the Fermi wave vector (k — kp). In particular, the be- 
havior of the self-energy and spectral function at k = kp 
determines the low-energy properties of the system. For 
k = kp there is only one peak in A{k, to) at u> = 0, which 
is a (5-function peak because Im£(fc.F, ui) oc aj 2 ln|w| as 
uj — > 9(. Thus, the doped (or gated) graphene is a 
Fermi liquid because it possesses a Fermi surface (or a 
discontinuity in momentum distribution function) whose 
presence is indicated by a 5 function in A(kp, uj) at uj = 
(or, equivalently, non-zero many-body renormalization 
factor). 

In Fig. 3 we show a quasiparticle lifetime due to 
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electron-electron interactions, which is given by l/r(k) = 
2Im[£(fc,£fc)]. Since the velocity of the quasiparticles 
close the the Dirac points is approximately a constant, 
the inelastic mean free path £ is obtained by £(£) — 
vft(£). In inset of Fig. 3, we provide the correspond- 
ing £, which shows that at n = 10 13 cm" 2 a hot elec- 
tron injected with an energy of 1 eV above Ep has an £ 
due to electron-electron interactions that is on the or- 
der of 20 nm. This will have implications for designing 
any hot electron transistor type graphcnc devices. Note 
that the calculate scattering time is a smooth function 
because both plasmon emission and interband processes 
are absent. In parabolic band semiconductors § plas- 
mon emission and interband collision thresholds cause 
discontinuities in the scattering time. 
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FIG. 3: Inelastic quasiparticle lifetime, 1/r = 2Im£(fc,ffc), 
in graphene due to dynamically screened electron-electron in- 
teractions, as a function of energy at T = for different den- 
sities. Inset shows the corresponding quasiparticle mean free 
path for n = 10 13 cm" 2 . 



SUMMARY 



In summary, we calculate the electron-electron interac- 
tion induced quasiparticle self-energy in graphene layer 
within the leading-order dynamical-screening approxima- 
tion. Our calculated spectral function indicates that 
the extrinsic graphene (i.e. Ep ^ 0) is a Fermi liq- 
uid. We also calculate hot carrier inelastic scattering 
due to electron-electron interactions. The linear disper- 
sion and chiral property of graphene give lifetime energy 
dependences that are qualitatively different from those of 
parabolic-band semiconductors. We can apply our calcu- 
lated self-energy to obtain other quasi-particle properties 
of graphene. 

This work is supported by US-ONR and LPS-NSA. 



[1] K. S. Novoselov et al, Science 306, 666 (2004); Nature 

438, 197 (2005). 
[2] For a review see A. K. Geim and K. S. Novoselov, Nature 

Materials 6, 183 (2007), and reference therein. 
[3] Y. Zhang et al, Nature 438, 201 (2005). 
[4] P. R. Wallace, Phys. Rev. 71, 622 (1947). 
[5] See, e.g., G. D. Mahan, Many Particle Phystcs, 3rd ed. 

(Kluwer/Plenum, New York, 2000). 
[6] E. H. Hwang and S. Das Sarma, Phys. Rev. B 75, 205418 

(2007). 

[7] E. H. Hwang, B. Y. K. Hu, and S. Das Sarma, 

|cond-mat/070 3499 (2007). 
[8] R. Jalabert and S. Das Sarma, Phys. Rev. B 40, 9723 
(1989). 

[9] S. Das Sarma, E. H. Hwang, and W. K. Tse, Phys. Rev. 
B 75, 121406 (2007). 



